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fvj . Abstract 

It is shown that during expanding phases of flat homogeneous cos- 
mologies all small enough non-linear perturbations decay exponentially. 
This result holds for a large class of perfect fluid equations of state, but 
O ■ notably not for very "stiff" fluids as the pure radiation case. 

cr 
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bJQ. 1 Introduction 

On expanding phases of nearly homogeneous cosmological models it is believed 
p\, • that small perturbations are pulled apart and so washed away. Thus, if they 

were the only type of perturbations present they would decay to zero, making 
the model more and more homogeneous as time elapses. This intuitive picture 
of homogenization, basic in all arguments used on inflationary models, has so 
far not been rigorously justified beyond the corresponding one on the linearized 
equations off a homogeneous background, the corresponding nonlinear case of 
Newtonian cosmological models [[[] , or the Electro- Yang-Mills- Vacuum case us- 
ing conformal methods, gj. 

The main mathematical difficulty to tackle this problem has been the absence 
of a set of variables in which Einstein's equations coupled to a fluid could at the 
same time have a well posed formulation and a linearization off a homogeneous 
expanding background where all eigenvalues have negative real part, for in that 
case, as we shall argue in §[| decay of small (but nonlinear) perturbations can 
be shown to follow quite easily. As we shall show in §g this problem has been 
recently overcome by using a novel set of Lagrangian variables, in which the time 
component of a frame variable follows the fluid four- vector under evolution, H. 
Thus, using these variables we are able to show exponential decay of small 
non-linear perturbations for a wide range of equation of state. 
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In the next section we introduce Friedrich's system for the Einstein-Euler 
system of equations in Lagrangian variables. Since this system has a well posed 
initial value problem (it is symmetric hyperbolic) we know local in time solutions 
exist, and so the only question is whether they decay or not if they are sufficiently 
close to an homogeneous expanding cosmology. To answer this we look at the 
linearization of this system off a homogeneous cosmological model, and display 
the conditions under which the real part of the eigenvalues of the system are 
negative definite. For the flat case, the only one considered in detail, they turn 
out to be only conditions on the equation of state, and the only one which seems 
to be a real restriction within the class of equations of state usually considered 
is condition ( p9[) , which limits the sound speed of the fluid to be less than one 
third of the speed of light. This rules out from consideration the stiff or pure 
radiation perfect fluids. As noted by several colleagues, perhaps this limiting 
case could be treated using the conformal Einstein's equations J2| . It seems that 
one can not go beyond that limit with the present approach and thus one is led 
to wonder whether this limitation has a physical origin, in which case it should 
be explored in more detail. 

In §H we complete the assertion on decay by giving a general proof of expo- 
nential decay for symmetric hyperbolic systems whose eigenvalues have negative 
definite real part, that is for system as the one we are considering. This theorem 
is a simple adaptation to the case under consideration of a more general class 
of results which are known in the literature on non-linear decay, for a resent 
review of them see || . 

We present the results for flat backgrounds, for non-flat ones similar results 
should hold as long as the expansion is large enough, but further work is needed 
to fully understand how big this expansion has to be in terms of invariant 
quantities, in particular how negative the cosmological constant has to be, so 
we don't elaborate in detail here. 



2 The evolution equations 

We use the evolution system for the Einstein-Euler equations as defined in j| . In 
that system one introduces a set of orthonormal frame vectors {u a , e"}, i = 1..3, 
u a u b g a f, = 1, efe^Poft = — %• In accordance with the Lagrangian coordinate 
scheme, the time-like component of the frame is chosen to be the four- velocity of 
the fluid, (so in general it is not surface orthogonal) , and the other frame vectors 
are chosen to be dragged along the integral curves of u a , that is along the flow 
lines, using Fermi transport, that is, To 2 *,- = 0. The remaining derivatives of the 
frame vectors are grouped into two set of scalars, a 1 :— IVo and Xij '■— ~^i°j, 
they shall form part of the evolution system. Note that Xij 1S n °t symmetric, 
for u a is in general not surface orthogonal. 

Part of the fluid equations are implied by the evolution equations for the 
components of the acceleration, a*, and a constraint which we shall indicate 
latter. The rest of them, and the only one surviving as such is the conservation 
equation for the energy density p := T a bU a u b , equation (||). 

To complete the evolution system one adds part of the Bianchi identities 
to obtain evolution equations for the components of the electric and magnetic 
decompositions of the Weyl tensor. The resulting system is, 



ftef = -Xi 3 ej ~ dieo (1) 

d t r{ k = -Ti kX i l -2 Xi[l a k] h jl -B ip e^ k (2) 

d t ai + v % D kX i k = -{xS+PxW+^tfi-xSVaj (3) 

dtXij ~ D j a i = ~Xi l Xlj - a i a j + 2ly2 X X[ij] ( 4 ) 

k p A 

dtp = -{p + P)X (5) 

d t E t3 + D k B l(l e 3) kl = X(i k E j)k + 2x k (iE j)k -h iJ x kl E kl -2 X E ij (6) 

K I 

+2aiE lk {l B j)k - -(p + p)(X(ij) ~ ^Xhij) 
d t Bij - DkEi^Ej) ' = x (i B ])k - xBij - XkiB pq e p ^e q ;) (7) 

~2aie (j-Ej-)fe 



where 



p + piTp _ _ .-, 
dp 2 



/ 3: =— ^773-^- ( 8 ) 



This system is symmetric hyperbolic when the speed on sound (y ) is pos- 
itive, and so in that case it has a well posed initial value formulation. Thus 
there is an energy estimate and local in time solutions exist and are unique for 
smooth enough initial data sets in L 2 . In what follows we shall use implicitly 
the existence and uniqueness of these solutions. 

For simplicity we are considering here only homentropic fluids. In the more 
general case one needs further evolution equation for the entropy and its space 
derivatives, see [|3|. The evolution equation for the entropy just asserts it is 
conserved along integral curves of the fluid four-vector, the evolution equation 
for the space derivatives of the entropy, assert consequently that they decay 
at the precise rate at which the expansion smoothes things out. Since these 
evolution equations are not coupled at the principal part level, one can always 



treat them as we shall treat the frame in §3.1 and obtain decay properties 
without having to impose further conditions on the system. 

In addition to these evolution equations there are a number of constraint 
equations, the usual ones that arise in general relativity, plus some connecting 
the different fields, for instance those which guarantee that the Xij are part 
of the connection coefficients of the frame. Here we give a brief discussion on 
the constraint problem, but do not go in details, for the aim of this work is to 
study the time behavior of solutions once initial data for them -satisfying the 
constraint- are given. 

To solve the constraints one would follow the usual procedure for frame equa- 
tions, and so we only remark on the differences which appear in this case. One 
first solve for a three metric and extrinsic curvature the usual scalar and vector 
constraints with some boundary or asymptotic conditions. Here we consider 
either an initial surface whose background is a flat 3-torus, or a flat R 3 . In 
the case of a 3-torus we consider periodic boundary conditions [J. In the open 

1 Here the perturbation of the mass density has to be adjusted so that there can be solutions. 



case (R 3 ) we consider isolated perturbations in the sense that one requires the 
perturbations to decay to zero asymptotically and are of finite energy. In the 
general case the perturbations in R 3 decay too slowly for the frame perturba- 
tions to be in L 2 , for there is a mass perturbation, and so one must use the 



variant of the decay theorem discussed in §3.1. One can also tune the mass 
density perturbation so that the total mass perturbation vanishes, but these 
seems to restrict the space of allowed perturbations. For this subclass of per- 
turbations one can use the standard decay theorems given in §|3[ Once this is 
done one then provides a lapse-shift pair which in this case can be thought of 
as the 3+1 decomposition of the fluid four vector. Thus one has a four met- 
ric, and one knows, using the field equations, all space and time derivatives of 
that metric at the initial surface. One then chooses a frame for that metric so 
that the time-like component of that frame coincides with the fluid four-vector. 
Most of the time derivatives of that frame are fixed from the metric condition 
and so are given by the extrinsic curvature of the surface, the rest, amounting 
to frame rotations and accelerations are fixed by the gauge condition imposed 
(Fermi transport) on the frame and from the acceleration, which is given by 
the Euler's equation, (jig), which here is part of the constraints. Thus we have 
all derivatives of the frame at the initial surface and so can compute the frame 
coefficients and from them all frame related quantities which enter as dynamical 
variables. They satisfy automatically the rest of the constraints. In the same 
way one can compute the Weyl tensor at the surface and its frame components, 
and so complete the set of initial data. 

One does not need for this case to prove that the constraints are propagated 
for this follows from general arguments, see for instance ||. Indeed suppose 
initial data is given satisfying the constraint, and suppose for contradiction that 
at some point p inside the domain of dependence of the initial slice the con- 
straints cease to be satisfied. Then we can take a small neighborhood of that 
point and inside it, to the past of it, construct a space-like slice so that p is in 
its domain of dependence. In that slice the smooth induced initial data satisfies 
the constraints and so we can generate a local solution of Einstein's equation 
using, say, the harmonic gauge. That solution satisfies the constraints along 
evolution. Furthermore we can construct locally a frame whose time-like com- 
ponent points along the fluid and whose other components are propagated using 
Fermi transport. The proof of this involves only the theory of ordinary differ- 
ential equations. Thus we have a solution of Einstein's equations in the gauge 
of our propagation equations and so it is a solution to them, furthermore all 
constraints are satisfied, even at p. But the solutions to the evolution equations 
are unique, and so we reach a contradiction, implying the constraints must prop- 
agate correctly inside the domain of dependence. The constraint propagation 
has also been checked explicitly, |I] 

In the above system it appears v 2 in the denominator, so it is not immediate 
that this equations behave nicely in the limit when this quantity -the sound 
speed- goes to zero, and for many equations of state this happens when the 
density goes to zero. Here we are not interested in reaching the limit where 
the expansion makes the density go to zero, but nevertheless it is interesting to 
obtain bounds which are uniform in v 2 in the whole region where it becomes 
small. An easy way to do this is to rescale the acceleration with a factor v and 
so get a system whose coefficients are bounded even in the limit p — > 0. Defining 
o-i — h a i, an d multiplying equation (0) by - we get, 



<9*ef = -Xi Je j - vcne^ (9) 

dtTi = -Tlx/-2v Xl[l & k] W l ~B v s™ k (10) 

d t a t + iyD kXl k = -(Xi j + /W +v 2 (xh-X l J ))a J (11) 

dtXij-vDjai = -x/xij - ("f + v 2 )aia 3 +2v 2 xXM ( 12 ) 

<9 f p = -(p + p)x (13) 

dtEij + DuB^erf 1 = X (i k E 3)k + 2x k (l E 3)k -h l3 x kl E kl -2xE t3 (14) 

K 1 

+2vai£ lk {l B 3)k - -{p + v){X(;i ) - -^Xhj) 
d t B l3 - D k E l(l e 3) kl = x k { l B 3)k ^ X B lJ -XkiB pq e p \ l e qk l) (15) 

-2vaie \iE 3 ) k 

where ^(p) :— — p ^Ji j-jr and J3 :— (3 — 7 = —37 — i^ 2 , to get these equations we 
have used the evolution equation for p, d5|) , and one of the constraint equations 
for the system, namely one of the fluid equations in the Euler description, 

D jP =(p+p)a 3 (16) 

This system is symmetric hyperbolic even for v — 0, and so well adapted 
for the case in which the mass density tends to zero. But to have a well posed 
system one also needs that the coefficient of the system be smooth. It can be 
seen that for many realistic equation of state 7 is bounded as a function of p, 
but in general it is not differentiable. For instance, for an equation of state of 
the form p — kp^ , 1 < 7 < 3 second derivatives are not bounded when p — > 0, 
unless the relation between p and p becomes linear, [j But notice that using the 
evolution equation for p one can see that time derivatives of 7 are bounded [for 
there is a contribution of an extra factor (p + p)]. Similarly, using (|l6|), one can 
see that the space derivatives of 7 are also bounded n. So, in general one can 
not proceed to make a priori estimates for derivatives of the fields, -as needed 
for showing existence- for equation (|lj) is known only a posteriori to hold. But 
nevertheless one can make a posteriori estimates on the solution, for if we know 
that the constraint are satisfied we can, using ([l6]) , obtain bounds on the space 
derivatives of the solutions to any desired order. This shall be enough for our 
purposes of obtaining uniform (in v 2 ) estimates in regions where solutions are 
known to exist. 

2.1 The background solution 

We take as background solution a Friedmann-Robertson- Walker Universe, namely, 

ds 2 = dt l-?°Vt( dx 2 + d y 2 + dz 2 ), (17) 



2 For an equation of state of the form p = kp^ one has, 7(p) = „ (1 + kp^ x ) 

3 For the case p = kp^ one can see that all derivatives are polynomial expressions in p 2 
and so regular, as long as 7 > 1. 



with u> = 1 + fcr 2 /4, and ao(t) solution of, 



,ao 2 k _ ftpo A 
a ag 3 3 



(-) 2 + 9-T + T = 0> (18) 



where fc = 1,-1,0 determines the curvature of the surface, A < is the cosmo- 
logical constant, and po is the background energy density. 
The background density satisfies, 

po = -3h (po+p ), (19) 

where ho — I 2 -, is the Hubble function and po — po(po) is the pressure, which 
is assumed to be a function of the density only. We shall assume that either 
A < or that we are in a time period where ho > ho > 0, for some constant ho- 
lt is easy to check that all evolution and gauge equations are satisfied if in 
the above coordinates one takes the following frame: 

o^ = o^ = (1,0,0,0), oe f = (0,-5f) (20) 

a 

The corresponding fields for this frame take the form: 



o^jk = ir^i-hjkx' + S'jXk), (21) 

za 

cTi j = — oXij = — ^o hij (22) 

all others = 0. (23) 

where x l = e t > 1/ = -^Stx". 

The proper time of co-moving observers would be taken to be the time with 
respect to which we assert the exponential decay result. Given any perturbed 
space-time this time is no longer an invariant, but it is so up to order e, that is up 
to the size of the perturbation. Furthermore, while the expansion is producing 
an exponential decay the difference between this time and the proper time of 
geodesic observers is uniformly bounded and goes to zero as e — > 0. 

2.2 The linearized equations 

The linearized equations are, 



ie" = -ho ief - \Xi J oe^ - voichSq (24) 

irj fc = -ho iFJfc - e l \ x B iX (25) 

iOj + v DjiXi 3 = -hoaiCLi (26) 

lXij ~ v a Djihi = -2h (iXij + '^o iX[ij]) - l-^ij (27) 

iBij — Dfe(iS TO ( i )£j) m = — 3h iBij (28) 

k 1 

tEij + D k { 1 B m(i )eo ]) km = -2>h iE tJ - - (p + p ){iXij - 3 iX%)( 29 ) 

pi = -3/i (l + ^o)/°i - (Po+Po) iX (3°) 



where 



= 1 - Zvl - 97, 

2.3 Negativity of the eigenvalues 

We notice that our system has the following form, 

u t = A a (u)L> a M+(B + eBi(M))M, (32) 

where the matrix A a {u) is symmetric, in the sense that there exists a strictly 
positive symmetric bilinear form H(u), smooth in u, such that H(u)(A a (u)) — 
(A a (u))^H(u) = 0. 

Due to this property it follows quite easily that the eigenvalues of the system 
have strictly negative real part for small enough e if and only if HBq + B* Q H 
is a negative definite bilinear form. In the next section wc shall show that the 
strict negativity of the eigenvalues suffices to show exponential decay. 

It can be easily seen that a family of bilinear forms making the system 
symmetric-hyperbolic is given by, 



< (e, T, a, x, E, B, p), H (e, T, a, x, E, B, p) > (33) 



C^hijh h I fem l ln — C y n didj 



1 



' x\- \X-ij Xmn "r XijXmn) ~r qX J 

+h im h jn E ij E mn + h im W n B l3 B mn + C 2 p p 2 

where xW (xtj) denote the symmetric trace free, and the antisymmetric parts of 
Xij respectively, Z M „ := ah^ — 2 u^qUv, and C e , Cr, C x , and C p are constants 
which can take any nonzero value, reflecting the fact that the bilinear form which 
symmetrizes the system is not unique. These constants shall be determined 
bellow in such a way as to maximize the range on which the negativity of Bq 
holds. Notice that H(u) — H$, a bilinear form which only depends on the 
background solution and so it is constant in space directions. 

Using the values of Bq from the linearized equations we obtain, 



-i- < (e, T, a, x , E, B, p), (H B + Bj5#b)(e, I\ a, x , E, B, p) > (34) 

-CJoLvhVefe" + C^h kl h mn r km P 



1 
a 



C^-afa, - 2C*[(1 - 3, 2 )( X A ) 2 + {x ST f + ^X 2 } 



-3(E) 2 - 3(B) 2 - 3C 2 (1 + v 2 )p 2 



~(%B lj e* i k hi P h im h kn r* n 

n 



;C2(l-3v*)+C 2 Jp + Po ))px 



-[C 2 + ^( Po +po)}E x ST 
i 

1 ,K 



V6 * v 

Thus negativity of the eigenvalues follows if we can choose the constants on the 
diagonal of H so that: 



1 C 2 



4a 2 C2 



< K (35) 



< h 2 (36) 



1 i2 „ 1,2 



3 "0- 9 n.,2 7^2 < ! ( a> °) ( 38 ) 



4^ < K (37) 

P o + Po d 2 p 
2v 2 dp 2 

Zvl 1 (39) 

1 [Cx+ <Pl±lA ] < , (40) 



2V6 L x 2C X 

r «(l + +3i/g)C x (po+po)C, 



2VT+^ 1 6C P C 



< fto (41) 



We now choose the constants in such a way to obtain the desired bounds. 
The first and the second can always be satisfied by choosing C e small enough. In 
fact, from the arguments used for the case where the frame perturbations are not 
square integrable one sees that these to conditions are superfluous. Note that 
the second is implied by the first if the fifth inequality is taken into account. The 
third can be always satisfied by choosing C'r small enough. The fourth inequality 
is a real constraint on the equation of state of the fluid, as is the fifth, which in 
particular rules out from our considerations pure radiation Universes, both are 
very similar. The sixth inequality can be maximized choosing C 2 — K,(po+po)/2, 
and the seventh, choosing C p = -4=. Both give similar conditions, the first 
gives, 

K(p +p )<12hl, (42) 

while the second, 

<Po + Po) < 12 ^"°} h 2 < 12h 2 (43) 

and so the second implies the first, and both are satisfied if (B9|) holds and 

K{p + Pa )<8h 2 (44) 



Using now (p_8|), and again ( |39| ) to estimate po < Po/3 in the second condition 
we get, 

k <(^-£). (45) 



a(t) 2 v 6 3 ; ' 



Since the right hand side is positive, while k = —1,0, 1, we see that this is 
only a condition for the case of closed cosmologies, namely k = 1. In that case 
it says, 

1 = k<1 2^-l A ^ (46) 

where Mq, and Rq are the mass and radius of the Universe. Thus slowly ex- 
panding Universes with mass over radius close to closure and small cosmological 
constant can not be treated with our method. 
We summarize this in the following Lemma: 

Lemma 1 // the equation of state is such that conditions (pq-ptf) are satisfied, 
and for the case k — 1 condition j[(q ) is also satisfied, then HqBq + BqHq is 
negative definite. 

3 The Energy Argument 

In this section we recall the standard argument leading to asymptotic stability 
of solutions from conditions on the non-principal part of a symmetric hyperbolic 
system. 

We begin considering a system of the form: 

u t = (A a + eAl{u))D a u + (B + eB 1 {u))u (47) 

with: 

1. The matrices Aq, and Bq are constant in space and time directions. 

2. Symmetric Hyperbolicity: there exists a strictly positive symmetric bilin- 
ear form H = H{u) = Hq + eHi(u) smooth in u, such that 

H(u)(A% + eA\{u)) - (>4g + eAt{u))^H{u) = 0. 

3. Decay condition: The matrix Bq is negative definite, namely 
< u, (HqBq + BqHoju >< —5 < u, H Q u >, for some S > 0. 

Remarks: 

1. Perturbations of flat Friedmann- Robertson- Walker space-times can almost 
be cast in this form by splitting the solution as a background solution plus 
a small part, ue '■= uo + su, and noticing that the background solution 
(out of which one constructs Aq, and Bo) is constant in space directions. 
In this case, the background solution, as well as Aq, Bq, and H depend 
on time, and so we shall latter modify the present argument accordingly. 

2. For our system H = Hq, for it can be chosen only to depend on the 
background fields, and on the tetrad components of the metric, which by 
definition are constant (in space) scalars. 

3. For our system, and only for convenience in notation, we shall later use 
derivatives along the frame directions and not D a . 



4. The topology of the initial space can be either compact, T 3 , say, or open 
i? 3 . In the presence of constraints the open case might be uninteresting 
for there could be too few solutions of the constraints which are in the 
required function spaces for showing existence. We shall make latter a 
variant of the standard theorems which requires less stringent conditions 
on the initial data. We shall always assume we have solutions of the 
constraint equations in the required functional spaces. 

We now define an energy vector by: 

E a := {d t ) a <u,Hu> + <u, HA a u >, (48) 

and compute, 

\7 a E a = <u u Hu> + < it, Hu t > + < u, H t u > (49) 

+ < D a u, HA a u > + < u, HA a D a u > + < u, D a {HA a )u > 

= 23? < u, H(u t + A a D a u) > + < u, (H t + D a (HA a ))u > 

= 23? < u, H B u > +e[23? < u, {H + eH^Bm > 

+ < u, {(HJt + D a {H Al + H x A a + eHxAiy)u >] 

Similar expressions hold for space derivatives of u, (i.e. tangent to the family 
of hypersurfaces) since due to the constancy along space and time directions of 
Aq and Bq, they satisfy similar equations with same zeroth order (in s) terms. 
Thus, considering a E" sum of all the E a, s for u and all its derivatives up to 
order p > 3, integrating in space and time its divergence, using Gauss theorem, 
and (for simplicity) taking the limit when the integration region shrinks to zero 
on the time direction we get, 

-(E p )<--E p + eF(E p ), (50) 

where 

E p (u) := / E a p n a {u) dS, (51) 

n a is the normal to a family of space-like hypersurfaces parametrized with r, and 
where we have used standard Sobolcv and Gagliardo-Nircmbcrg-Moscr estimates 
to control all the nonlinear terms. 

Thus, given an initial value for u such that its E p norm is finite we can 
choose e so that the right hand side of the above inequality is negative, and so 
we conclude that the E p norm of u at future times can not become larger that 
its initial value. Furthermore if the function F(-) can be chosen to approach 
zero at least linearly then for small enough e there is an exponential decay. But 
the above mentioned condition in F(-) follows directly from the smoothness 
assumptions on the coefficients of the equations system, thus we have: 

Theorem 1 Given any initial data uq whose E p norm is finite, then there exists 
£o > such that for all < e < Eq a solution exists and it decays exponentially 
to zero. 
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For the case at hand one must be a bit more careful for the frame chosen 
is not surface orthogonal, and furthermore the matrices Aq are not time inde- 
pendent. To control the first problem we shall use the background homogene- 
ity hypersurfaces for which the zeroth component of the background frame is 
surface-orthogonal, thus the contributions due to the non-surface-orthogonality 
of the frame to the energy difference will be of order e. We call that time r 
(= to), and correspondingly choose n a := oe°. We also introduce the matrix A 1 
so that the principal part becomes 

u t = {A + eA\(u))e t (u) + {B + eB x (u))u (52) 

For the case at hand, namely the flat case, this change is trivial, and in particular 

Bo = -Bo- 
Defining now the energy four- vector as 

E a := eg < u, Hu > +ef < u, BA l u >:= e^E° + e a L E\ (53) 

we get a similar expression for eqn. (E9I) , 

V a E a = <u t ,Bu> + <u, Bu t > + <u, B t u > (54) 

+ < u, HA i e i (u) > + < u, e^HA^u > + < e,(u),Mii > 

+ E°W a e a Q + PVaef 

= 23ft < u, B(u t + A l e t (u)) > + < u, (B t + e^BA^u > 

+ E° X + E\r%-a t ) 

= 23ft < u , B B u > +E°ox + e[23ft < u, (B Q + sB x )B x u > 

+ < u, ((#i) t + ei(H A\ + H x Ai + eH x A\))u > -E 1 ^ - l0i )] 

where we have used that o^i = 0, and for the flat case o^ ki = 0. The only new 
zeroth order term is E° x = 3h Q E°. Since, d^a^^gE^) = a^ 3 y^gV a E a - 
ao 3 v /= 5 oX E °, w e get 



d^V^E") (55) 



.-3 / — - "P 



a , .. 

= < u t .Bu > + < u, Bu t > + < u, B t u > + < u, BA l e l (u) > 
+ < u^e^BA^u > + < e i {u),HA i u > +E i \7 a e°; 
23ft < u, B{u t + A l e % (u)) > + <u,(B t + e,(BA l ))u > 
+E i (ri ji -a i ) 

23ft < u, BoBou > +e[23ft < u, {H + eB x )B x u > 
+ <u, {{H x ) t + e,(B A\ + H x Al + eH x A\))u > -#(il^ - l0i )] 

Thus applying Gauss theorem to the coordinate divergence in the above 
expression, and noticing that a^ \J—g = 1 + 0(e) we get the desired energy 
estimated as for the four dimensional flat case treated above. Notice that the 
new volume element is time independent to zeroth order, and so the constant 
on the Sobolev embedding and on other estimates can be taken to be also time 
independent. 
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Theorem 2 If the equation of state of the fluid is such that conditions ( p& 
p{\ ) are satisfied, then given any initial data perturbation of an expanding flat 
homogeneous cosmology, uq whose E p norm is finite, there exists So > such 
that for all < e < Eo a solution exists and it decays exponentially to zero. 

3.1 A variant for the case in which ef is not in L 2 . 

In some cases it is not possible to have perturbations satisfying the constraint 
equations and having the L 2 norm of ef finite, this is because in the open 
(R 3 ) fiat background case there could arise the need of considering massive 
perturbations, and they do not decay at infinity sufficiently fast. In that case 
we can only have, D a e^ £ L 2 . There are several ways to deal with this problem, 
and perhaps the best treatment is the one in |7|. For the present problem 
a machinery such as in M is not necessary, for the decay properties can be 
established locally (pointwise) for part of the equations. The idea is not to look 
at the Sobolev norms of e^ 1 , but rather directly at its maximum norm of it, this 
is possible because the time derivative of the frame does not have any space 
derivative, it is just a ordinary differential equation. Thus, if wc call by e the 
frame variables, and by u all the others we have a system of the form, 

e t = -h e t + fu{u)+ef e (e,u) (56) 

Ut = {Al+EA a 1 {e,u))D a u+{B Q + EB l {u))u 

where here the matrices A a , and B are not the same as before, but are the corre- 
sponding restrictions. The usual theory of energy bounds using the Gagliardo- 
Niremberg-Moser estimates, (see for instance Proposition 3.7, and note bellow 
it, of 0) give us for the second equation, and the corresponding equations for 
the space derivatives of both, e, and u, the following inequality, 

4-(E p (u) + £ P -i(£>e)) < (57) 

CLT 

-6{E p (u) + S p _i(£>e))+eF(£; p (u),S p _i(£)e),||e||co), (58) 

for some smooth F with F(0, ■, •) = 0, since the negativity of the restricted Bq 
as well as for the whole B (needed for the estimates on the system of space 
derivatives (De,Du)) follows as before. Note that the ||De||co is bounded by 
the Ep^i(De) norm if p > 3. 

While for the first equation we get, at each point q of the base manifold, 

|e(?,*)|t < -S\e(q,t)\+G u (E p (u))+EG e (E p (u),\e(q,t)\), (59) 

With G u , and G e smooth, polynomially bounded, functions of all its arguments, 
and with G u (0) — 0, and G e (-, 0) = 0. In fact in our case G u is a linear function 
of the E p {u) norm. 

To show exponential decay for e small enough we assume, for contradiction, 
that given initial data (eo, uq) at t = there exists a time T* > given by 

C E p (u(T)) = E p (u(0))(l + A 2 ), or 

T* = mi : I S p _i(Ue(T)) - £ p _i(£>e(0))(l + A 2 ), or ), (60) 

T> ° I INIco m = ||e|| C o (0) + ^M)} 
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where C(E p (u)) is a function to be defined bellow, and A is some non-zero 
constant. We shall show that T* = oo, and from the proof it also follows that 
the solution not only exists for all times, but in fact decays exponentially. 
Given the initial data we choose sg small enough so that 

e G F(E p (u(0))(l + A 2 ),E p ^(De(0))(l + A 2 ), ||e|| c o (0) + g^MM ) < ( 6 1) 

6 -(E p (u(0)) + E p ^(De(0))) 

then it is clear that (u, De) will decay as e~~ as long as T < T* 
We now turn to the frame equation and choose El so that 



s L G e (E p (u(0))(l + A 2 ), ||e|| c o ( o) + C ^ m ) < (62) 

6 C(E p (u(0))) 
2 CI l e llc°(o) H g )) 

then we have, as long as T < T* 

\e(q,t)\ t < - S -\e(q,t)\+G u (E p (u)), (63) 

and from this it follows that 

IMIco(T) < (l|e||c»(0) + j) e ~^> ( 64 ) 

where C — G u (E p (u)). Thus we see that taking C — C, and s = mm{sQ,EL} 
we obtain that T* = oo, and so the solution exists for all times and furthermore 
decays exponentially. Thus we have proven a similar Theorem as the one above 
for the case where the L 2 norm of the frame is not bounded, but rather its C° 
norm is. 

Acknowledgements: 

I thank H. Friedrich for pointing to me the possibility of using his descrip- 
tion of the Einstein-Euler system for this particular problem, and for several 
enlightening discussions. I also thank Uwe Brauer and Gabriel Nagy for point- 
ing to aspects of earlier versions of this work which needed further refinements 
or corrections. 



References 

[1] U. Brauer, A. Rcndal, and O. A. Reula, "The cosmic no- hair the- 
orem and the nonlinear stability of homogeneous Newtonian cosmo- 
logical models," Class. Quant. Grav. 11 (1994), no. 9, 2283-2296, 
http://xxx.lanl.gov/gr-qc/9403050 

[2] H. Friedrich "On the Global Existence and the Asymptotic Behavior of So- 
lutions to the Einstein-Maxwell- Yang-Mills Equations," J. Diff. Geometry 
34 (1991), 275-345. 



13 



[3] H. Friedrich "On the evolution equations for gravitating ideal fluid bodies 
in general relativity," Phys. Rev. D 57 (1998), 2317-. 

[4] H. Friedrich "Personal Communication" . 

[5] Reula, O.A. "Hyperbolic Methods for Einstein's Equations," Living Re- 



3reula/| ), 1998 



Relativity, http://www.livingreviews.org/Articles/Volurnel/1998- 



[6] Kreiss, H-0, and Lorenz, Jens "Stability for time-dependent differential 
equations," Acta Numerica (1998), 203-285. 

[7] Choquet-Bruhat, Y., Christodoulou, D. and Francaviglia, M. "Cauchy data 
on a manifold," Ann. Inst. Henri Poincare XXIX, no. 3, (1978) 241-255. 

[8] Taylor, M. E. "Partial Differential Equations III" Applied Mathematical 
Sciences 117, Springer, 1997. 



11 



